Three-dimensional aspects of fluid flows in channels. II. Effects of Meniscus and Thin 

Film regimes on Viscous Fingers 
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We perform a three-dimensional study of steady state viscous fingers that develop in linear chan- 
nels. By means of a three-dimensional Lattice-Boltzmann scheme that mimics the full macroscopic 
equations of motion of the fluid momentum and order parameter, we study the effect of the thick- 
ness of the channel in two cases. First, for total displacement of the fluids in the channel thickness 
direction, we find that the steady state finger is effectively two-dimensional and that previous two- 
dimensional results can be recovered by taking into account the effect of a curved meniscus across 
the channel thickness as a contribution to surface stresses. Secondly, when a thin film develops in 
the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement 
with experimental results. The effect of the thin film renders the problem three-dimensional and 
results deviate from the two-dimensional prediction. 
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I. INTRODUCTION 

Interfacial instabilities in three-dimensional channels 
give rise to a rich phenomenology in systems that range 
from nano and microscales[l[ to macrometric channels [2, 
[H, 0] , and from which a number of practical applications 
can be drawn. 

For instance, controlled drop breakup in micro- 
channels has proved useful in the fabrication of low poly- 
dispersity micro-emulsions |5l and in the enhancement 
of micro-reaction processes [g, Q. In the latter, three- 
dimensional effects are crucial, as they are responsible 
of a vortex flow structure within the droplet Q that en- 
hances the mixing process of the reactants. 

A widely studied interfacial instability in channels is 
that of fingering, which occurs whenever a low-viscosity 
(or high-density) fluid drives a high- viscosity (or low- 
density) one. The instability, first studied by Saffman 
and Taylor 0, leads to interface dynamics where finger- 
like structures emerge and compete. The problem has a 
steady-state solution, composed by a single finger of con- 
stant velocity U and occupies a fraction A of the width 
of the channel. 

Experimentally, finger growth has been studied mainly 
in Hele-Shaw cells. These consist of a pair of plates of 
length L and width W separated by a thickness b. For 
such systems, it has been pointed outfioj that the station- 
ary finger is determined by a single control parameter, a 
modified capillary number defined as l/B = 12Ca/e 2 . 
For a fluid with viscosity r\ and surface tension a, the 
capillary number, Ca = rjU/a, measures the competi- 
tion between driving forces, such as viscous stresses and 
gravity, and restoring forces, like surface tension. l/B 
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also includes the degree of asymmetry of the cell, given 
by the aspect ratio e = b/W. If l/B is the only con- 
trol parameter of the system, all experimental data, i.e. 
all finger widths, should be described by a single curve 
when plotted as a function of this parameter. Contrary 
to this view, experiments show that there exists a fa mily 
of curves X vs. l/B for different aspect ratiosfll], Il2{ . 
This fact suggests that a three-dimensional effect, given 
by the interplay between the dynamics in the channel- 
thickness and in the channel- width, is determinant for 
the steady-state solution. 

Theoretically, fluid-flow in a channel at small veloci- 
ties pertains to the lubrication regime, in which the flow 
occurs mainly along the direction of L given that it is 
much larger than both W and b. Hele-Shaw flows are 
a limiting case in lubrication theory, where b is much 
smaller than W. Owing to the smallness of b, the prob- 
lem is rendered effectively two-dimensional by averaging 
all fields over the thickness of the channel. Averaging the 
equations of motion also reduces the interface from a sur- 
face to a line, often called the leading interface. In views 
of the averaged model, three-dimensional effects enter as 
perturbative corrections to the boundary conditions that 
hold at the leading interface in terms of Ca and e, par- 
ticularly to the Gibbs-Thomson condition, which relates 
the pressure drop across the interface to the interface 
curvature and surface tension. 

Progress towards a three-dimensional description of 
the problem has been made since the pioneering work 
of Saffman and Taylor Q, who solved the problem of a 
stationary finger in the absence of surface tension in two 
dimensions. McLean and Saffman[l3| included the ef- 
fect of surface tension and were the first to obtain a 
A vs. l/B prediction by solving numerically the two- 
dimensional model. According to their results, A is a 
monotonically decreasing function of l/B that saturates 
to A — > 1/2 as l/B — > oo. The prediction of McLean and 
Saffman is unique in l/B, so the role of the aspect ratio 
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is precluded from their theory. 

The relevance of three-dimensional effects was first sug- 
gested by Park and Homsv[T3|. who pointed out that a 
thin film of fluid in the channel-thickness direction would 
contribute to the pressure drop at the leading interface. 
Using perturbation methods for slightly curved leading 
intcrfaces(small e), they found that for low Ca the pres- 
sure drop varies as Oa 2 / 3 ,a result that matched the early 
prediction of Bretherton[15j for capillary tubes. 

Sarkar and Jasnow[l6j used the modified pressure drop 
to solve the steady state finger. Their results agreed bet- 
ter with experiments but were restricted to low values 
of Ca. It was shown by Tabeling and LibchaberpH 
that corrections to the pressure drop can be used to 
reduce three-dimensional experimental data to the two- 
dimensional results of McLean and Saffman. A modified 
pressure drop can be accounted for as an effective sur- 
face tension. Using the correction of Park and Homsy, 
Tabeling and Libchabcr were able to reduce their data 
to McLean and Saffman results for moderately low val- 
ues of 1/-B, where fitting parameters were used to esti- 
mate the correction terms. In an experimental study [l2T|. 
Tabeling, Zocchi and Libchaber observed that, contrary 
to the McLean-Saffman prediction, the finger width can 
go below the one-half limit for sufficiently high 1/ B and 
sufficiently large e. 

Reinelt extended the expansion of the pressure drop 
up to O(l) in Ca and included the effect of larger as- 
pect ratiosp"?]]- Computation of the steady state finger 
yielded solutions that better agreed with experiments for 
0(1) values of Ca. For small e Reinelt observed a better 
agreement between numerics and experiments. However, 
for relatively large e this agreement is lost. 

Higher Ca values have only been explored in the case 
of flat leading interfaces by Halpern and Gaver[l8j . Their 
numerical results are consistent with results found by 
Reinelt and Saffman (l9j for Ca 0(1) and e = 0, and 
show that the pressure drop is insensitive to the capil- 
lary number for Ca > 20. 

As an alternative to the sharp interface model, a num- 
ber of mesoscopic approaches have gained importance 
in interface dynamics. These are based on order pa- 
rameter evolution equations of the Cahn-Hilliard type. 
Being mesoscopic in nature, fluids are separated by dif- 
fuse regions instead of sharp interfaces, where the inter- 
face boundary conditions arise naturally. All mesoscopic 
models that address the viscous fingering problem so far 
are two-dimensional. For fluids of arbitrary viscosities 
and densities, Folch et a/[2(| [2l[ used a set of coupled 
evolution equations for the velocity potential and order 
parameter that describes accurately the early stages of 
destabilization of the leading interface, and approaches 
McLean and Saffman results as the viscosity of the dis- 
placing fluid is made negligible. The strict one-sided sit- 
uation, were one of the fluids is inviscid, was studied by 
Hernandez-Machado et al^^. They used a single evo- 
lution equation for the concentration that includes dy- 
namic effects in the form of chemical potential gradients 



and described correctly the steady state finger. 

In a preceding paper (23[, we have shown that a detailed 
three-dimensional description of fluid-flow in a channel 
can be done by means of a mesoscopic model which 
we implement numerically via a Lattice-Boltzmann al- 
gorithm. The model considers a fluid-fluid interface in 
contact with solid boundaries. In contrast to classic ap- 
proaches, it allows for slip at the contact line by means 
of a diffusive mechanism inherent to the mesoscopic na- 
ture of the interface. This circumvents the complica- 
tions of contact line dynamics in the classic formulation, 
associated to the viscous dissipation singularity [24j. In 
Ref . [23l] . we focused on the case of a flat leading inter- 
face. We showed that depending on the velocity of the 
contact lines, which we control by modifying the diffusion 
strength, the interface can either advance as a meniscus 
or develop as a finger. In the latter case, a thin film of 
displaced fluid is left adhered to the walls of the channel. 

In this paper we extend our Lattice-Boltzmann simu- 
lations to the case of a non-flat leading interface, where 
a viscous finger is expected to appear. Our aim is to 
provide a detailed description of the mechanisms that 
affect the steady finger and that cause deviations from 
two-dimensional results. To do so, we study fingers that 
form in the meniscus and thin film regimes separately. 
We cover values of Ca up to O(10) and explore various 
aspect ratios. 

The paper is organized in the following manner. In 
Scc.|IT]wc present the governing equations of the system 
which we solve numerically by means of a Lattice Boltz- 
mann algorithm, presented in the preceding paper [2^]. 
Results are presented in Sec. IIIII In Sec. IIII Al we de- 
scribe the simulation strategy and parameter steering 
procedure. As a validation test, in Sec. IIIIBI we com- 
pute the dispersion relation of the interface in the two- 
dimensional limit and compare it to the analytic predic- 
tion of the Saffman- Taylor problem. Sec. IIII Cl is devoted 
to the study of stationary viscous fingers; in Scc lIII C Tl wc 
focus on fingers pertaining to the meniscus regime in the 
channel thickness, which we found to be effectively two- 
dimensional, while in Sec. IIII C 21 fingers in the thin film 
regime are studied. We find that fingers in the thin-film 
regime are three-dimensional and cannot be described 
by the two-dimensional theory in general. A discussion 
of our results where we compare with previous experi- 
ments is presented in Sec lIVI In Sec. [Vj we present the 
conclusions of this work. 



II. GOVERNING EQUATIONS 

We consider the motion of two viscous fluids, whose 
dynamics are governed by the Navier-Stokes equations, 

p (d t v + v{V ■ v)j = —VP - (jNpL + riV 2 v + pg. (1) 

Here v is the fluid velocity, P is the pressure, p is the 
density, r\ is the fluid viscosity and g is the acceleration 
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FIG. 1: Schematic representation of the system. Dashed lines 
indicate projections of the fluid-fluid interface in the xy and 
xz planes. The leading interface corresponds to the xy pro- 
jection. 



of gravity. The extra term, 0V/i, is mesoscopic and ac- 
counts for intcrfacial forces. Here, <j>(r, t) is an order pa- 
rameter and p(4>) is the chemical potential. <f> has the 
property of being uniform in the volume of each phase 
and non-uniform in an intcrfacial region of typical size £. 
In the present case, volume values are chosen as <j> = ±1 
for the displacing and the displaced fluid respectively, so 
the interface is located at <j> = 0. The size of the interface 
is set to £ = 0.57. 

The dynamics of (j> obey a convection-diffusion equa- 
tion, 



d t <t> + v-V<l> = MV>, 



(2) 



where M is a mobility coefficient. In equilibrium, the 
pressure and chemical potential minimize a free energy 
functional, from which explicit expressions P(p) and p((f>) 
can be derived. For further details the reader is referred 
to Ref.[H]. 

We work on a linear channel, composed by two solid 
plates of width W and length L parallel to the xy plane, 
separated by a distance b, as depicted in Fig. [T] There 
exist two principal directions in the system: a lateral di- 
rection, parallel to the xy plane, and a transverse one 
parallel to the xz plane. We will denote these by sub- 
scripts || and _!_ respectively. 

The impenetrability and stick boundary conditions 
at the walls are enforced by setting v{x, y,z = 0) = 
v{x,y,z = b) = and <pv(x,y,z = 0) = <f>v(x,y,z — 
b) = 0. At both ends of the channel in the x direction 
the flow is homogeneous. Hence, d x pv(x = 0, y, z) = 
d x pv(x = L, y, z) = 0, and d x <j)v(x = 0, y, z) = d x (j)v(x = 
L,y,z) = 0. Periodic boundary conditions are imposed 
in the y direction. 

As for the fluid-fluid boundary, the Gibbs-Thomson 
relation is recovered by integrating Eq.|T]) across the in- 
tcrfacial region, 
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where a is the surface tension and R a is the radius of 
curvature of the interface in the direction a. 

We now briefly review the classic treatment of the 
problem. First, v± is assumed to be much smaller than 
v\\, which in turn is expected to be parabolic in z. As a 



result, Eq.lU) is recast in the form of an average veloc- 
ity field which holds in the volume of each fluid, called 
Darcy's Law, 
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where triangular brackets denote an average over the 
channel thickness. Under these conditions, R± is ex- 
pected to be much larger than Ru . Hence, in the two- 
dimensional theory the Gibbs-Thomson relation is sim- 
plified to AP = a/R\\. 

Corrections to this expression arise whenever 1/R± 
is not neg ligible. For such cases, Libchaber and 
TabclingpJl have proposed that thin film effects can be 
accounted for by defining an effective surface tension 



a* = a 1 
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so the two-dimensional form of the Gibbs-Thomson con- 
dition is recovered. For this purpose, they used the esti- 
mation of Park and Homsy[lJ| of the pressure drop for 
Ca — > and slightly curved leading interfaces (e — > 0), 
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(5) 



As a result, their experimental results collapsed to the 
McLean- Saffman curve when using the corresponding 
definition of the control parameter, 1/B* = (a/a*)l/B. 

We solve numerically Eqs. {1} and ([2]) by means of a 
Lattice-Boltzmann algorithm. For further details of the 
method, the reader is referred to the preceding paper [2^]. 



III. RESULTS 
A. Simulation Parameters and Setup 

The traditional description of the viscous fingering 
problem corresponds to situations in which the relevant 
forces at play are viscous stresses and capillarity. For 
the particular case of fingering in a Hclc-Shaw cell these 
forces are expressed in terms of a modified capillary 
numberpjl 1/B = \2{W/b) 2 {AriU + A pgb 2 / 12) / a, where 
A?7 and Ap are the differences in viscosity and density 
between the fluids. 

To ensure that capillarity and viscous forces dominate 
the dynamics of the fluids, inertia must be small com- 
pared to both of these forces. Wc enforce this situation 
by neglecting the convective term in Eq. ([1]) . As for com- 
pressibility, we consider low Mach number flows, which 
we achieve by keeping U <C c s . For our scheme, it suffices 
to set U < 0.01. 

Our goal is to explore the viscous fingering problem for 
a wide range in l/B. Due to computation resource limi- 
tations, e is restricted to O(10) at most for the majority 
of runs. To achieve high values of 1/B, say O(10 3 ), Ca 



4 




0.5 1 1.5 2 2.5 

k' 

FIG. 2: Dispersion relation for the linear stability of the 
interface. Simulation parameters (in simulation units) are 
a = 0.046, g = 0.1, b = 11.0 for all runs; (+) Apg = 3.3 x 10~ 6 
and (■) Apg = 6.6 x 10" 6 

must then be O(10). Our strategy is to keep the interface 
velocity and the viscosity in ranges of U ~ O(10~ 2 ) and 
r/ = O(10 _1 ). Hence, Ca can be tuned by means of the 
surface tension. 

The channel is implemented as follows. We set a 
rectangular simulation box of dimensions Nx x Ny x 
Nz. Due to the flow symmetry, we simulate only one 
fourth of the real channel by setting boundary condi- 
tions as follows: d y pv y (x,y = 0, z) = d y pv y (x,y = 
W/2, z) = 0, d y 4>v y (x,y = 0,z) = d y 4>v y (x,y = 
W/2,z)0, d z pv z (x,y,z = 0) = d z <j>v z (x,y, z = 0) = 0. 

B. Linear Stability in the two-dimensional limit 

We first verify the linear stability of the interface, 
i.e., the behavior of an initially flat interface that has 
been subjected to a small perturbation. We study fluids 
of equal viscosities, so the instability is gravitationally 
driven. This is done by fixing the body force term in 
Eq. |T]) as pg = 8r]/b 2 U exp ((j> + l)/2, where U exp is the 
maximum expected velocity for a Poiscuillc flow. The 
modified capillary number reduces to 1/B = W 2 Apg / 'a, 
In this case, the linear stability analysis of the interface 
evolution of the averaged equations yields the dispersion 
relation 

= ^KApg ~ ^ 2 ), (6) 

where u> is the exponential growth rate of a sinu- 
soidal perturbation to the flat interface solution. The 
perturbation is characterized by its wavelength, A = 
lix/k. By considering dimensionless frequencies cu' = 
uj/(U/2W)B 1 / 2 and wavenumbers, k' = WB x l 2 k, the 
dispersion relation becomes universal, i.e., u>'(k') = 
k'(l-k' 2 ). 

We prepare a base flow corresponding to a flat inter- 
face in the xy plane that propagates at constant velocity. 
The interface in the xz plane is nearly flat throughout the 



simulation, so the system is effectively two-dimensional. 
Once the base flow is fully developed, the interface is 
shifted according to a single mode perturbation of wave- 
length A = W and an initial small amplitude. We follow 
the evolution of the amplitude, A(t), which is measured 
as A(t) = Xti p (t) —x(t), where xu p and x(t) are the inter- 
face tip and mean interface positions respectively. The 
growth rate, 10, is extracted as a linear fit of log(A(t)) vs 
t. 

Fig. [2] shows a comparison between the universal dis- 
persion relation and our results. To quantify the degree 
of accuracy of these results, we fit our data to the gen- 
eral form ak'(b — ck' 2 ). We find a most unstable mode 
at k' max ~ 0.56 and a first unstable mode at k' ~ 0.96, 
both 4% below the exact result. 



C. Viscous Fingers 

In a preceding study p3j. we have shown that it is pos- 
sible to control the generation of a thin film in the channel 
by adjusting the diffusivity of the order parameter. Al- 
though for usual experimental conditions this is not a rel- 
evant parameter (it might be relevant in nano-channels), 
it gives the possibility of elucidating the role of a thin 
film in viscous fingers. Diffusivity is accounted for by a 
Peclet number, Pe = Ub/D, where D is the diffusion co- 
efficient. By combining the effects of Pe and Ca, one can 
either suppress or induce the formation of a thin film. In 
particular, a small value of the product CaPe results in 
the suppression of thin films, while the contrary is ob- 
tained for high CaPe. Results from the preceding work 
give a penetration threshold of CaPe ~ 10 _1 . 

The strategy is to study first fingers for which CaPe < 
10 _1 and then extend this simulations to CaPe ^> 10 _1 . 



1. Meniscus Regime 

We first study fingers for which no film of displaced 
fluid develops in the xz plane of the channel. We carry 
out simulations with modified capillary numbers in the 



TABLE I: Control parameters and finger width for runs in of 
the meniscus regime. 



Run 


e 


Ca 


CaPe 


1/B 


A 


(a) 


0.17 


0.11 


0.08 


99 


0.709 


(b) 


0.17 


0.22 


0.16 


198 


0.675 


(c) 


0.17 


0.45 


0.19 


290 


0.640 


(d) 


0.06 


0.11 


0.04 


522 


0.558 


(e) 


0.06 


0.19 


0.02 


1045 


0.525 


(f) 


0.06 


0.23 


0.03 


1672 


0.523 


(g) 


0.06 


0.48 


0.11 


2090 


0.529 


00 


0.06 


0.68 


0.21 


3136 


0.518 


(i) 


0.04 


0.74 


0.26 


4175 


0.521 


(j) 


0.05 


0.76 


0.27 


6012 


0.519 
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FIG. 3: Interface snapshots at two different times for e = 0.17(the plot is off-scale), 1/B — 99 and CaPe = 0.08. The thick 
line parallel to the xy plane corresponds to the leading interface, while the thick line parallel to the xz plane corresponds to 
the interface projection in the center of the channel. Thin lines correspond to the contact lines. The first snapshot corresponds 
to t = O.llb/U, while the second snapshot, at t = 17.746/{/, corresponds to the steady state finger. 
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FIG. 4: Collapsed interface profiles in the xy plane for the 
meniscus regime. Parameter values corresponding to each 
symbol can be consulted in Table [T] The error (small bar at 
the right bottom) corresponds to one lattice spacing. The 
larger bar indicates the size of the diffuse interface, approxi- 
mately 3£. 



range 100 < 1/B < 6000. We have studied different ge- 
ometries, ranging from e = 0.17 to e = 0.04. The aspect 
ratio is decreased by decreasing the channel thickness. 
We summarize the simulation parameters in Table [U 

For each run we observe the usual phenomenology for 
the leading interface. During the early stages of interface 
evolution, the amplitude of the perturbation grows until 
a finger emerges and widens. This stage is followed by a 
relaxation of the interface shape, until a Saffman- Taylor 
finger develops. The finger propagates with a steady ve- 
locity U, leaving behind a growing region where the finger 
has flat sides. In this region a constant finger width, AW ', 
can be defined. As for the channel thickness, we observe 
that the initially flat interface rapidly relaxes to a menis- 
cus, which also has a steady shape. In Fig. ([3]) we show 
a three-dimensional plot of the interface for run (a) in 
Table [J at two different times. In the plot we show both 
the contact lines and the leading interface; both contact 



lines follow the leading interface. 

To check for consistency in the steady state solution 
we use the semiempirical interface profile obtained by 
Pitts (25[, which reproduces experimental results accu- 
rately for a wide range of finger widths. The equation 
for the interface shape reads, 

cos(7n//2A) = cxp(7nz;'/2A). (7) 

where x' and y' measure the distance from the finger tip 
in units of half the channel width. The natural scalings 
in this equation are nx' /2A and ny' /2A. Consequently, 
all profiles should collapse into a single curve if these 
scalings are used. Fig. U shows plot of interface profiles 
corresponding to runs of Table [IJ As expected, all inter- 
face profiles fall in the same universal curve within error. 
In addition, our collapse is in fair agreement with Eq.([7]). 

The selection rule in the viscous fingering problem 
is expressed as the functional dependence of the finger 
width with the modified control parameter. We compare 
our results with the A vs. 1/B prediction of McLean 
and Saffman. We find that runs with e = 0.17 show 
wider fingers than predicted, while runs with smaller e 
agree better with the two-dimensional result. Even in 
the absence of a thin film, the xz interface projection 
has a certain curvature. This can be accounted for by 
defining an effective surface tension in terms of the radii 
of curvature of the interface, which we are able to mea- 
sure directly. The effective surface tension then reads 
a* = er(l + Rn/R±). The correction factor in this ex- 
pression is given by the quantity in parentheses, which 
increases for strongly curved meniscus. The rescaled con- 
trol parameter then reads 1/B* = (1/B)/(1 + R\\/R±). 
Of course this correction should be more evident in the 
low 1/B region, where A varies rapidly with the modi- 
fied control parameter. In Fig. [5] we show a plot of A vs. 
1/B*. Points fall on the McLean-Saffman curve for the 
wide range of 1/B* considered, regardless of the aspect 
ratio. 
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FIG. 5: Finger width as a function of the rescaled control 
parameter 1/B* in the meniscus regime. 
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FIG. 6: Interface projections in the xy and xz planes for runs 
with different CaPe values. Plots correspond to the same 
simulation time, (a) CaPe = 0.85, (b) CaPe = 4.44. 



2. Thin Film Regime 

We now extend our simulations to fingers in the film 
regime. Penetration in the xz plane occurs for high 
CaPe, so we choose to sample 1/B at fixed D. Conse- 



TABLE II: Control parameters and finger width for thin film 
regime runs. 



Run 


r 


Ca 


CaPe 


1/B 


A 


(a) 


0.25 


2.80 


12.32 


835 


0.592 


(b) 


0.25 


3.36 


17.74 


1002 


0.589 


(c) 


0.25 


6.61 


68.61 


2004 


0.569 


(d) 


0.25 


15.9 


400.41 


4003 


0.558 


(e) 


0.35 


8.96 


1515 


1403 


0.549 


(f) 


0.49 


34.7 


4330 


5247 


0.527 


(g) 


0.64 


50.9 


3973 


5247 


0.517 


(h) 


0.78 


68.5 


7192 


5247 


0.508 


(i) 


1.00 


91.9 


8019 


5247 


0.493 


0) 


1.00 


131 


156598 


5430 


0.494 



qucntly, CaPe increases with increasing 1/B. To resolve 
the thin film correctly we must take into account the 
finite size of the interface. As explained in Ref.[23j], the 
thin film is insensitive to the channel thickness already for 
b = 23. We therefore choose sufficiently thick channels. 
We explore a wide range of aspect ratios, 0.25 < e < 1.0 
and modified capillary numbers, 800 < 1/B < 5300. 

We first explore the CaPe ~ 0{1) region, close to 
the penetration threshold. In Fig. [6] we show interface 
projections in the xy and xz planes located at z = b/2 
and y = W/2 respectively. We show two sets of in- 
terfaces, corresponding to two different CaPe values; 
(&)CaPe = 0.85 and (b)CaPe = 4.44. In Fig.^a) the in- 
terface in the xz plane presents a penetrating structure, 
but a well developed film is absent. The finger in the 
xy plane is not well developed either, and it presents an 
anomalous tip. Conversely, in Fig. [6jb) both interface 
projections describe well developed fingers. It is then 
clear that deviations from the Saffman- Taylor finger in 
the xy plane are correlated to the interface structure in 
the xz plane. An interesting feature of the run corre- 
sponding to Fig.[6ja) is that the the xz interface structure 
is persistent. This means that the length of the finger in 
the xz plane is constant in time, a consequence of the slip 
velocity of the contact line. The diffusion strength is not 
large enough to maintain a meniscus, which in the one 
hand makes the slip velocity smaller than the channel ve- 
locity. Nevertheless, as the interface relaxes to a thin film 
shape, curvature deviations from equilibrium increase the 
slip velocity, making the contact line advance to restore 
the meniscus shape. 

We next explore the range CaPe > O(10) for which 
simulation parameters and observed finger widths are 
summarized in Table ILT1 In Fig. [7] we present snapshots of 
the three dimensional interface at two different times for 
run (b) in Table HT1 The first snapshot corresponds to the 
early stage of the finger formation. Looking at the inter- 
face projections in the xy plane, we see that the contact 
linc(light line) is close to the leading interface(dark line) 
and no film is present in the xz plane. In the next snap- 
shot the contact line has moved away from the tip, thus 
giving rise to the growth of a wetting film. The shape 
of the finger is in agreement with the typical morphol- 
ogy found in experiments. To illustrate this, in Fig. [5] we 
compare the shape of the finger to Eq. ([7]). Within error, 
our profiles are consistent with Pitts shape. 

Fig. [9] shows the measured finger width as a function 
of 1/B. The lowest aspect ratio we have considered cor- 
responds to e = 0.25(runs (a)-(d) in Table ITT]) . We see 
that for all 1/B values considered the finger width falls 
above the McLean-Saffman prediction. We increase the 
aspect ratio to e = 0.35(run (e) in Table Hlj) . As a result, 
the measured finger width decreases. Runs for which e is 
larger confirm this tendency in a systematic way. Tests 
(f)-(j) in the same table correspond to a fixed value of 
1 /B with increasing e. We find that for sufficiently large e 
the finger width goes below the one-half theoretical limit 
of McLean and Saffman. 



FIG. 7: Interface snapshots at two different times for e = 0.25(the plot is off-scale), 1/B — 1002 and CaPe = 17.74. Thick 
lines correspond to the xy and xz interface projections in the center of the channel. Thin lines correspond to the contact lines. 
Times are t = 0.57b/U and t = 28.846/(7. 
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FIG. 8: Rescaled interface profiles for the thin film regime. 
Symbols correspond to data presented in Table HIl The bars 
in the bottom at the right indicate the error bar and diffuse 
interface size as in Fig. [4] 
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FIG. 10: Finger width as a function of the rescaled control 
parameter for the thin film regime. 



IV. DISCUSSION 

Our results show that the finger width decreases with 
increasing aspect ratio. To maintain 1/B fixed while 
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FIG. 9: Finger width as a function of 1/B. 



varying the aspect ratio of the channel, one has to vary 
Ca accordingly. As a consequence, the film thickness and 
the capillary pressure are altered. If we increase the as- 
pect ratio(as in the high-l/i? region in Fig. then Ca 
must decrease to keep 1/B fixed. As a consequence, the 
film thickness and the capillary pressure decrease as e 
increases, which is consistent with a narrower finger. 

This behavior has been observed, for instance, in ex- 
periments by Tabcling, Zocchi and Libchabe r|l2ll . and 
addressed in numerical calculations by Rcinelt[l7l] where 
the effect of the thin film was introduced perturbatively 
in the two-dimensional model. Experiments suggest that, 
for high 1/B, increasing the cell aspect ratio has a thin- 
ning effect on the finger, which is what we observe in 
our simulations. Results of Rcinclt suggest the opposite 
tendency. 

The aforementioned experiments were done at small 
e and Ca, and at high viscosity contrast, defined as 
c = (r]2 — i]i)/{r)2 As we have shown in Ref . (23| . 

the thin film thickens as c — ► 1. Under these conditions, 
experiments show that the finger width goes below the 
one-half limit even for cells with e = 0.009. This is due 
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to the small thickness of the film which is a consequence 
of the low Ca and high c values used in experiments. 
In our simulations the thin film is about t/b ~ 0.25, 
much thicker than the experimental value, t/b ~ 0.05. 
As a consequence, curvature effects in our simulations 
are strong enough to keep the finger width above one 
half even for high values of e. To achieve the experimen- 
tal regime thinner film should be considered. We have 
considered a cell aspect ratio of e ~ 0.05 and c = 0.9. 
Nevertheless, Ca is still too large, the film is then thick 
enough to keep us in the low 1/B* regime, where the 
finger width is still larger than one half of the channel 
width. Due to computational limitations we do not ex- 
plore smaller e. 

The fact that for a given 1/B there exist different fin- 
ger widths for different aspect ratios raises the doubt of 
l/B as being the only control parameter present in the 
system. To this end, we compute the rescaled surface 
tension a* = a (1 + Rn/Rx), where the radii of curva- 
ture are measured at the finger tip. We then rescale the 
control parameter according to 1/B* = (a/a*)l/B. In 
Fig. |TD] we show a plot of the finger width as a function 
of the rescaled control parameter. At low 1/B*, results 
agree with McLcan-Saffman results. We conclude that in 
this region the finger is effectively two dimensional and 
that three-dimensional effects can be accounted for even 
at Ca ~ 1. 

At high values of the rescaled control parameter, our 
results deviate systematically from the McLean- Saffman 
curve, until the finger width goes below the one half 
limit of the two-dimensional theory. This behavior is 
qualitatively different from the one found for the menis- 
cus regime, in which the McLean-Saffman curve could 
be recovered at any value of 1/B* . Hence, we conclude 
that deviations from two-dimensionality are caused by 
the thin film. 

An important feature in the A vs. 1/B* plot is that 
finger width appears to be determined by 1/B* uniquely. 
This suggests that 1/B* is the only control parameter of 
the problem. 

We have explored a region of values of the aspect ra- 
tio between the Saffman- Taylor(e — > 0) and Rayleigh- 
Taylor(e = 1) limits of the fingering instability. In both 
limits, the relevant control parameter appears to be an 
effective modified capillary number. In addition, the in- 
terface shape is remarkably universal, as suggested by 
FigsHandE] 

V. CONCLUSIONS 

We have carried out a detailed study of the viscous 
fingering problem in three-dimensional channels for fluids 
of different densities and viscosities. We have studied the 



single finger solution for systems in which either a thin 
film develops across the channel thickness or a meniscus 
is stabilized. 

For systems in which no thin film is present, McLean- 
Saffman two-dimensional results describe the dependency 
of the finger width as a function of a rescaled modified 
capillary number, 1/B*, which has a correction that de- 
pends curvature of the interface direction of the channel 
thickness. This holds for arbitrary high values of 1/B*, 
evidencing that a complete displacement across the chan- 
nel thickness renders the problem two-dimensional. 

We have extended our studies to situations where a 
thin film develops across the channel. We find different 
values of the finger width when changing the channel as- 
pect ratios at fixed modified capillary number, an obser- 
vation that is consistent with previous experiments [l2| . 
This non-uniqueness seems to disappear as the control 
parameter is corrected by curvature effects associated to 
the thin film, i.e., when the finger width is compared to 
1/B*. 

For low 1/B*, the finger width collapses to the 
McLean-Saffman curve. However, at high 1/B* the fin- 
ger width deviates from this curve, and goes bellow the 
limit of one half in units of the channel width. 

Our work is done at high values of the capillary num- 
ber. Consequently, the effective capillary pressure in 
our simulations is large enough to keep the finger above 
the one-half limit of the two-dimensional theory for high 
values of the channel aspect ratio. Experiments in 
Refs.fTD. [HJ were done in cells with e ~ 0.03 and at 
Ca ~ 10~ 3 , a regime that falls beyond the scope of this 
work for computational reasons. Nonetheless, for low 
1/B*, we recover the same results, indicating that the 
same mechanisms hold, even if the actual aspect ratio 
and capillary number are not the same in experiments 
and simulations. 

To our knowledge, experiments of fingering in high 
aspect ratio channels have not been conducted so far. 
Our results could be confirmed, for instance, in micro- 
channels, where the aspect ratio is typically large and 
in which the meniscus to thin film transition could be 
observed. This is then an open question for experimen- 
talists to confirm. 
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